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The flux of neutrino-induced muons at the surface of the Earth is calculated from injection of neutralino 
annihilation products in the core of the Sun and the Earth. An improved treatment of neutrino propagation 
through the Sun is performed and the results are presented in an easy-to-use parameterization. For an explicit 
supersymmetric model, an observable neutralino annihilation signal is demonstrated. 



1. Introduction 

The dark matter in the Universe may be con- 
stituted by the neutralino, which probably is 
the lightest supersymmetric particle. If these 
particles exist they will get trapped gravita- 
tionally by the Earth and the Sun where they 
can annihilate and produce other particles, e.g. 
fermion— antifermion pairs, gauge bosons and 
Higgs bosons. These particles can hadronise 
and/or decay producing high energy muon neu- 
trinos which can be detected in neutrino detec- 
tors, in which the neutrinos produce muons via 
charged current interactions and the muons can 
be detected due to the Cerenkov light they emit. 

We have considered the whole chain of pro- 
cesses from the annihilation products in the core 
of the Sun or the Earth to detectable muons at 
the surface of the Earth, see § for details. The 
neutrino flux at the surface of the Earth has been 
calculated earlier by others , e.g. but in this 

calculation the neutrino propagation through the 
Sun has been considered more carefully and the 
neutrino interactions near the detector have been 
simulated to get the muon fluxes. 

The hadronisation and/or decay of the annihi- 
lation products have been simulated with Jetset 
7.3 H and the neutrino interactions on the way 
out of the Sun have been simulated with Pythia 
5.6 M. The neutrino interactions near the detec- 
tor (in which the detectable muons are produced) 
have also been simulated with Pythia 5.6. 
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2. Annihilation channels and interactions 
of annihilation products 

A pair of neutralinos can annihilate to produce 
a fermion-antifermion pair or gauge bosons, Higgs 
bosons and gluons, i.e. xx ~ * ^ + ^~, qq, gg, qqg, 
W + W~, Z°Z°, Z°H°, W ± H^, H°H°. These 
annihilation products will hadronise and/or de- 
cay, eventually producing high energy muon neu- 
trinos. Since the density in the core of the Sun 
and the Earth is high, the possibility of interac- 
tions of the annihilation products with the sur- 
rounding medium must be considered. 

In agreement with earlier work [g] , we find that 
light hadrons and muons get stopped well before 
they have time to decay both in the Sun and in 
the Earth. Tau leptons lose hardly any energy at 
all and their interactions can to a good approxi- 
mation be neglected. Gauge bosons, Higgs bosons 
and top quarks also decay long before interacting. 
The gluon channels do not give rise to any signif- 
icant high energy muon neutrino fluxes, and can 
thus be neglected. The branching ratio to neutri- 
nos directly is close to zero for slow neutralinos 
annihilating and these channels do thus not con- 
tribute. However, the heavy hadrons (B's, D's, 
Ab's etc) may or may not interact and their in- 
teractions can not be neglected and they can not 
be considered to be stopped completely either. 

For the heavy hadrons we simulate the hadro- 
nisation as if the surrounding medium were not 
present. Afterwards we estimate how many times 
the heavy hadrons have interacted and how much 
energy they have lost in each interaction. The 
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Figure 1. a) The neutrino spectrum (y^ + v^) at the surface of the Sun from injection of 500 GeV 
neutrinos and anti-neutrinos in the core of the Sun with the method of Ritz and Seckel and with 
the new method, b) The neutrino spectrum [y^ + v^) weighted by E% at the surface of the Sun for 
the W + W~ -channel with a neutralino mass of 500 GeV. The two peaks in the Ritz and Seckel spectra 
correspond to and respectivley. 



neutrino is assumed to have lost the same frac- 
tion of its energy as the heavy hadrons. These 
heavy hadron interactions are taken care of both 
in the Sun and in the Earth. 

The Higgs bosons decay to ordinary particles 
and the muon neutrino flux they produce can be 
calculated from the flux that their decay products 
produce (l) . For the top quark, only the standard 
model decay t — > W + b has been considered. In 
supersymmetry, the decay mode t — > H + b is open 
if the H + mass is light enough, and the flux from 
this decay mode can be calculated from the fluxes 
of the b quark and the decay products of the H + . 

3. Neutrino interactions 

The produced muon neutrinos can interact on 
their way out of the Sun or the Earth. The total 
neutrino-nucleon cross section is given approxi- 
mately by ace = a{E v /GeV) x 10~ 39 cm 2 and 
&nc = b{E u /GeV) x 10 -39 cm 2 for charged and 
neutral currents respectively. We have calculated 
the coefficients a and b using the new GRV struc- 
ture functions || and a Q 2 -cnt of 0.3 GeV 2 . The 
result is a vn = 8.81, a up = 4.51, a p „ = 2.50, 
a Pp = 3.99, b vn = 2.20, b up = 1.97, h n = 1.15 
and byp = 1.14. 

One can easily find that neutrino interactions 
can be neglected in the Earth, but not in the Sun. 



The effective thickness of protons and neutrons, 
T p and T n , of the Sun is calculated by using the 
solar model of Bahcall et al. @ with the result 
T p = 1.1 xlO 12 g/cm 2 andT n = 3.6 xlO 11 g/cm 2 . 

We have approximated the Sun with a piece 
of homogeneous material with the thicknesses of 
protons and neutrons given above. For each pro- 
duced neutrino we have calculated the mean free 
path for interactions and if an interaction has 
taken place, it is chosen to be a neutral or a 
charged current. In case of a charged current, 
the neutrino is treated as absorbed and in case of 
a neutral current, the scattering is simulated with 
Pythia 5.6 and the procedure is continued with 
the new energy of the neutrino until the neutrino 
has reached the surface of the Sun. 

This differs from the approach of Ritz and 
Seckel [0 , which is the common approach, where 
the energy loss is considered to be continuous and 
neutral currents are assumed to be much weaker 
than charged currents. Neither of these assump- 
tions is very good. In general, a neutrino only 
participates in a few interactions (~0-2) on the 
way out of the Sun. Hence the process is highly 
discrete. In Fig. ^a, the neutrino spectrum at the 
surface of the Sun is shown for injection of 500 
GeV muon neutrinos and anti-neutrinos in the 
core of the Sun. The difference between the ap- 
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Parameterization (according to Eq. (fy)) of the muon fluxes for different annihilation channels. With these 
values the unit of the flux in Eq. (]l|) is m _2 (annihilation) _1 . An angular cut of 9 < 5° has been applied. 
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Annihilation in the Earth 
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proach of Ritz and Seckel and the new method is 
substantial. The two peaks in the Ritz and Seckel 
spectrum correspond to and respectively. 
In Fig. [l]b, the neutrino spectrum weighted by 
E^, is shown for a typical annihilation channel. 
The two peaks for and can still be seen. 
Ef,dN / dE u is proportional to the differential rate 
of events since the cross section for muon produc- 
tion (charged current) is proportional to the en- 
ergy and the range of the muon is proportional to 
the energy. The difference between the two meth- 
ods can clearly be seen, even though it has been 
somewhat washed out when applying the meth- 
ods to a neutrino spectrum in the core of the Sun. 
However, the difference in the high energy tail of 
the spectrum is significant and it gets important 
due to the ^-enhancement. If one uses the Ritz 
and Seckel approach, the error in the total rate 
is significant, and of the order of 5-20% too low 
with the higher error at higher neutralino masses 
(~ 1500 GeV). 

4. Simulation results 

The muon flux at the surface of the Earth has 
been calculated for different annihilation channels 
and for different neutralino masses. The annihi- 



lation channels considered are cc, bb, tt, t + t~ , 
W + W~ and Z°Z° for the reasons stated ear- 
lier. The simulation has been performed with 
neutralino masses of 50, 100, 150, 200, 250, 350, 
500, 750, 1000 and 1500 GeV and a top quark 
mass of m t = 150 GeV. For easy use, the calcu- 
lated fluxes have been parameterized as 
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where is the neutralino mass in GeV, z = 
E/j,/ m x an d Pit • ■ ■ )f>8 are the parameters fitted 
to the simulation results. The values of the fitted 
parameters are given in Table [l]. If one calculates 
the total flux above a certain threshold, these pa- 
rameterizations are accurate to ~ 15% as long as 
the threshold is not too close to the neutralino 
mass (Ejj 1 < 0.2m^). The parameterizations are 
poorest for the cc-channel and the 66-channel in 
the Sun at high neutralino masses (> 500 GeV), 
but these channels are not expected to dominate 
the flux at these high masses. The accuracy of 
the t + t~, W + W~ and Z°Z° parameterizations 
are usually of the order of 5% and never worse 
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Figure 2. The expected event rate, T, in con- 
centric angular bins (the first one corresponding 
to ±5° around the Sun and the next ones being 
of the same size in solid angle) for the supersym- 
metric model described in the text. The solid line 
corresponds to the total event rate and the dotted 
line is the background [ [l0[ . The errors shown are 
just the square root of the number of events in 
each bin. A cut of E'jf 1 > 10 GeV is also applied. 



than 10% and these channels usually dominate 
the flux. 

The event rate at a detector is given by 
^events = T ' aA c j f J2i B i®i where T A is the an- 
nihilation rate, A e ff is the effective area of the 
detector in the direction of the Sun/Earth, Bi is 
the branching ratio for annihilation channel i and 
$j is the muon flux (over threshold for the detec- 
tor) per unit area and annihilation. T A and B t 
depend on the supersymmetric parameters cho- 
sen and are considered in other papers (e.g. M). 

5. Example 

For a specific supersymmetric model and a spe- 
cific detector, one can calculate the event rate. 
Consider a neutrino detector with an effective 
area of 2000 m 2 (a typical size for Amanda || 
in the first set-up) in the direction of the Sun. 
For a Minimal Supersymmetric Model with pa- 
rameters n = 300 GeV, Mi = 600 GeV and 
tan j3 = 2.0, giving ~ 250 GeV, one gets the 
expected event rate shown in Fig. || where the an- 
nihilation rate and the branching ratios are given 
by || . Note that this specific example is not yet 
excluded by other experiments. 



6. Conclusions 

The treatment of neutrino propagation through 
the Sun has been improved (the approach of Ritz 
and Seckel Q was shown to give an error of up to 
~ 20%) and the calculated muon fluxes are given 
in an easy-to-use parameterization. The large 
neutrino telescopes now being built should have 
good possibilities to search through new parts of 
the supersymmetric parameter space. 
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